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I show that the usual model of the rotational response of a neutron star, which predicts rotation- 
induced neutronic vortices and no rot at ion- induced protonic vortices, does not hold (i) beyond a 
certain threshold of entrainment interaction strength nor (ii) in case of nonzero E~ hyperon gap. I 
demonstrate that in both these cases the rotational response involves creation of phase windings in 
electrically charged condensate. Lattices of bound states of vortices which are caused these effects 
can (for a range of parameters) strongly reduce the interaction between rotation-induced vortices 
with magnetic-field carrying superconducting components. 



Microscopic calculations of properties of interior of a 
neutron star suggest a presence of a mixture of superfluid 
neutrons, superconducting protons and normal electrons 
(see e.g. pj -0|). Where protons form a type-II super- 
conductor. This model suggests that the magnetic field 
in the star should induce a large number of vortices in 
the protonic condensate [l], 0] , while rotation should in- 
duce neutronic vortex lattice. Because the magnetic field 
configuration in a young star is very complicated the pro- 
tonic vortex structure is expected to be a complicated 
vortex tangle [3, 0] . It is important that in this standard 
picture for a typical star the average distance between 
protonic vortices is much larger than the magnetic field 
penetration length (which is of order of 10-100fm). Such 
a vortex tangle is indeed a thermodynamically unstable 
state but estimates suggest that a very large time scale 
is needed to expel these vortices 

The microscopic calculations also indicate that in a 
neutron star the effective mass of proton is very different 
from the bare mass @ . This implies that there is a strong 
dissipationless drag effect jjj: i-e. superfluid velocity of 
neutronic condensate drags superfluid density of protonic 
condensate and vise versa. In particular a neutronic cir- 
culation drags along some density of protons resulting in 
magnetic field generated by neutronic vortex. Therefore 
a neutronic vortex interacts magnetically with protonic 
vortices 0|- 

In the usual model for neutron star, the free energy of 
the mixture of neutronic and protonic condnensates has 
the form: 



Ppn 



where B = V x A is magnetic field, while v„ = 
(l/2m„)V(/>„ and v p = (1 / 2m p )V 4> p — (e/m p )A are the 
superfluid velocities of neutron and proton condensates 
in units h = c = 1 (generalizations to cases with non-s- 
wave pairing is straightforward) . Here m n ss m p — m are 
the bare masses of a neutron and a proton and <fi n ,p are 
the phases of the corresponding condensates. The third 
term in ([1]) represents current-current interaction Q . Be- 
cause of it the particle current of one of the condensates 



(w PlB ) is carried by the superfluid velocity of another: 

w p = Ppp v p + PpnVn] W„ = p nn V n + p pn V p , (2) 

where p pn — p np are the drag coefficients which can 
be expressed via effective masses m* n as follows p pn = 

Pnn m " m m " ■ Different microscopic calcula- 

~ 0.3m p to 
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tions 
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give for m* the values ranging 



0.9m p . Thus the drag strength can be as high as 
Ppn ~ —0-7p pp - From the eq.fl]) it follows that the elec- 
tric current induced by protonic circulation or by neu- 
tronic drag is given by: 
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When there is a 2-7T phase winding in cj> p (i.e. for an 
integral around the vortex core §V(j) p = 2vr), but no 
phase windings in 4> n (we denote this vortex as (1,0)) the 
circulation of electric charge is 
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In the constant density approximation the magnetic field 
obeys the London equation [§]: B + A 2 curl curlB = 
<I>o<5(r) where $o is the magnetic flux quantum and 
A = m/e^/ppp is the magnetic field penetration length. 
Imposing a 2ir winding in (f> p in eq. (|TJ) one obtains a 
vortex with finite energy per unit length: 
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where a is the cutoff length associated with vortex core. 

In the case when a 27r winding is imposed only on 
the neutronic condensates [lets denote it (0,1)], the drag 
effect produces some protonic current: 
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Indeed this drag current in turn produces magnetic flux 
which is not quantized but is a linear function of the drag 
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strength: ^ 



dlA — ££2l2jl (here the integration is 

done over the contour a located where J « 0) . The main 
contribution to the energy of (0,1) vortex is associated 
with the kinetic energy of the superflow of neutrons. This 
contribution is logarithmically divergent as a function of 
the system size R: 
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The secondary contribution is associated with the drag- 
induced kinetic energy of the protonic supercurrent and 
the magnetic energy. Because electrically charged cur- 
rents are localized at the length scale A (in contrast to 
neutronic superflow), this contribution is much smaller: 
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If a protonic fluxtube intersects with a neutronic vortex 
which carries co-directed flux, the magnetic energy and 
kinetic energy of supercurrents rise. The energy of the 
intersection can be estimated as: 
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According to [3|, |4| E l ch w 5M eV. Because protonic vor- 
tex configuration is complicated and because protonic 
vortices are much more numerous than neutronic vor- 
tices there should be strong pinning between the 
protonic and neutronic vortex matter. Much interest in 
the physics of vortex interaction was sparked by the re- 
cent calculations in [3| which suggest that the picture of 
strongly pinned protonic and neutronic vortices may be 
highly inconsistent with the slow precession observed in 
a few isolated pulsars. Although most of pulsars do not 
have such precession, this inconsistency casted doubts on 
the validity of the usual neutron star model. According 
to 0] , the precession would be possible only if the drag 
strength was orders of magnitude lower than values ob- 
tained in microscopic calculations. Following the work 
j3l| some alternative scenarios of magnetic response of a 
neutron star interior were proposed. In 12j, based on 
the calculations of £~ hyperon gap (l3| . it was proposed 
that one has a mixture of two oppositely charged conden- 
sates which gives further alternatives to type-I/type-II 
dichotomy [14j (for other aspects of the theory of two- 
component charged mixture see fl5L [lil . [l7j). The al- 
ternative conjecture put forward was that protons form 
type-I superconductor [3|, [9|, [Uj , suggesting that the mag- 
netic field would penetrate star interior via macroscopic 
normal domains. However it was pointed out in [Til ] that 
the energy difference of a neutronic vortex located in- 
side superconducting region versus inside a large type- 
I domain is quite large, which again leads to a strong 
interaction between neutronic vortices and flux-carrying 



domains. Let me stress that importantly the magnetic 
response of a type-I superconductor is highly nonuniver- 
sal. Generically this magnetic response involves forma- 
tion of multi-domains. The precise shape of the domains 
depends on a number of factors and is magnetic-history 
dependent [is| . Importantly a boundary between a nor- 
mal metal and a superconductor in magnetic field implies 
that there is a magnetic field and supercurrents within 
the range of A at this boundary. Thus even if protons 
form a type-I superconductor, neutronic vortices will 
be strongly interacting with the boundaries of normal 
domains which typically are plentiful in a type-I super- 
conductor [lij . The strength of the resulting pinning 
of neutronic vortices will depend on the domain's struc- 
ture (which is usually very complicated). Besides that 
the mechanism responsible for formation of large normal 
domains in a type-I superconductor is assocaited with 
the dominance of attractive core-core interaction between 
type-I vortices. However this force has the range of the 
coherence length. For a dilute vortex system with low 
mobility (as expected to be the case in a neutron star) 
the coalescence tendency of type-I vortices (and there- 
fore the tendency to form large normal domains) would 
therefore be extremely small. So in general, the interac- 
tion strength between neutronic vortex lattice and flux- 
carrying structures in type-I case should not be expected 
to be dramatically lower than the interaction between 
neutronic vortices and protonic vortices in a type-II case. 

Therefore it is an interesting question if there exists 
a mechanism which, in principle, could lead to any sig- 
nificant coupling reduction between neutronic and pro- 
tonic vortices. Here I propose two scenarios which lead 
to such an effect. First let me observe that beyond a cer- 
tain threshold of the drag strength, the usual neutronic 
vortices are not thermodynamically stable. That is, the 
vortex lattices form because such a state minimizes the 
free energy in a rotating system. The free energy of a 
(0,1) vortex in a system rotating with the angular veloc- 
ity CI is [| 
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Where M is the vortex momentum 
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I d\[- Ppn \{V<p n - e^A)| + - Pnn \V<p n \] (8) 
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In eq. ([8]) the first contribution comes from protonic cir- 
culation, while the second term corresponds to the neu- 
tronic circulation. The supercurrent of protons is con- 
centrated around the vortex within the range of pene- 
tration length (of order of 10 or lOOfm in the standard 
picture), while the characteristic size of the star is of 
order of 10 km. Therefore the protonic contribution 
/ d 3 r(r/m)p pn \V4> n - e(p pp / p pn )A)\ in the vortex mo- 
mentum (the last term in Eq. ([?])) is much smaller than 
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the neutronic contribution (which is not exponentially 
localized). Note also that this small contribution has 
the role of the energy penalty because p pn < and thus 
protons are circulating in the direction opposite to the 
neutronic circulation. Let us now rewrite the equation 
© as 
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Here we separated the energy associated with the elec- 
trically neutral superflow (given by the first term). It 
is associated with the phase gradient decoupled from 
A. The remaining terms describe the kinetic energy 
of the electrically charged currents (represented by the 
phase gradients coupled to vector potential A) and the 
magnetic energy. In a free energy of (0,1) vortex in a 
rotating frame, there is (i) a free energy penalty com- 
ing from the negative-drag-induced momentum of pro- 
tons, and (ii) energy penalties from the kinetic energy 
of protonic currents and magnetic field. Let me observe 
that the simplest vortex with phase winding only in the 
neutronic phase is not thermodynamically stable when 
\Ppn\l Ppp > 1/2. Thermodynamic stability is achieved 
when the free energy is minimized with respect to all the 
degrees of freedom which include phase windings. Previ- 
ously it was assumed that rotation induces vortices with 
only neutronic phase winding [2j. However from the eq. 
(|9|) we can see that if there is a 2n winding in </>„, and 
—Ppn/ Ppp > 1/2, the system can minimize the free energy 
in a rotating frame by creating an additional 2n winding 
in the protonic phase 4> p (around the same core). Lets 
denote such composite vortices (1,1). Its free energy is 
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Note that microscopic calculations [6[ estimate the drag 
strength in a neutron star can be as high as p pn / Ppp ~ 
—0.7. Let me observe that with increased drag strength 
beyond the threshold —p pn /Ppp > 1/2, the magnetic 
energy and the energy of protonic currents around a 
rotation-induced (1,1) vortex decreases. Moreover the 
additional phase winding does not affect the dominant 
neutronic negative contribution (coming from — MQ) to 
the free energy jJJ . Also because the protonic current 
of (1,1) vortex 
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has a circulation opposite to that of a (0,1) vortex 

j(o,i) K (JlAv^-eA 
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it carries an opposite momentum. Thus instead of a pos- 
itive free energy penalty — M • SI it produces a negative 
free energy gain in this term. Therefore Fr 1 ' 1 ^ < Fr 
for —p pn / Ppp > 1/2 and thus the (0,1) neutronic vortices 
are not thermodynamically stable in this regime. Instead 
a superfiuid supports rotation by creating a lattice of com- 
posite vortices (1,1). Observe that in general, for a given 
superfiuid momentum, there is also a possibility to min- 
imize the energy of charged currents for smaller ratios of 
\Ppn\l Ppp by creating a 2irN winding in <p n and compen- 
sating the contribution of V<p n in the second term in (J9]) 
by a 2ttM "counter- winding" in <f> p . But in the case of a 
vortex with a multiple winding in <p n the first term in ([9]) 
depends quadratically on TV while the last term depends 
on it linearly. Because in a neutron star, the p nn is very 
large compared to p pn and p pp the composite vortices 
with multiple windings in cf> n , like e.g. (2,1) should not 
be thermodynamically stable. 

Another physical consequence of the possible forma- 
tion of composite (1,1) vortices is that for —p pn /Ppp > 
1/2 with increased drag strength \p pn \ there is a decrease 
in the protonic current and magnetic energy contribu- 
tions to the energy cost of an intersection between the 
rotation induced vortex and protonic (1,0) flux tube. For 
nearly parallel vortices it is given by 
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Also the (1,1) vortices should have reduced interaction 
with the boundaries of the normal domains (if pro- 
tonic condensate is type-I). To estimate more accurately 
vortex-vortex interaction potential, it is important to ob- 
serve that the composite (1,1) vortex has a core in pro- 
tonic condensate due to protonic phase winding. This 
gives rise to the attractive core-core interaction with the 
range of coherence length [n| (which originates from 
minimization of the energy cost of suppression of the pro- 
tonic order parameters in the cores) . Therefore if protons 
form type-II condensate then for p pn / p pp > 0.5 the inter- 
action between approximately co-directed (1,1) and (1,0) 
vortices is further minimized by core-core interaction be- 
cause both (1,1) and (1,0) vortices have zeros of protonic 
condensate. 

There is another scenario leading to appearance of 
composite vortices and possible electromagnetic coupling 
reduction in a neutron star. It has been discussed that 
a star can have another electrically charged condensate 
associated with E~ Cooper pairs 12, Tj^. Let me observe 
that from the theory of two-component charged mixture 
[l^ | it follows that the presence of the second charged 
component alters qualitatively the rotational response of 
the system (if intercomponent Josephson coupling is neg- 
ligible |1G|). That is, in the free energy functional of such 
a system, one can form a gauge-invariant combination of 
the phase gradients cx (V(A, + ^s)) which is decoupled 
from the vector potential [lja, [l6| (here is the phase of 
hyperon condensate). It means that if another charged 
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component is present, then there is an additional super- 
fluid mode (i.e. besides the neutronic) associated with 
co-flow of oppositely charged neutronic and hyperonic 
Cooper pairs. Such co- flows transfer mass without charge 
transfer. From the calculations in 15, 16| it follows the 
system has two kinds of vortex excitations which carry 
identical momenta associated with these co-flows. These 
vortices have phase windings (A# p = 2tt, = 0) and 
(A9 P = 0, A#£ = 2tt). On the other hand these vortices 
carry different fractions of magnetic flux quanta: 
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where l^^x;)! an d n r i( P ,'s) are protonic and hyperonic 
condensate's number densities and masses. The density 
v of rotation-induced vortices in a mixture of c harg ed 
condensates depends on the angular velocty fl as [l(| 



(nip + m-s)Cl 



(15) 



Therefore these vortices should be more numerous than 
neutronic vortices (observe that this effect takes place 
no matter how small the superfluid density of the sec- 
ond charged component is). At the same time the 

phase stiffness of the composite superfluid mode is 
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deed much smaller than the phase stiffness of the neu- 
tronic condensate and therefore ordering energies of the 
rotation-induced protonic/hyperonic vortex lattice are 
much smaller than those of neutronic vortex lattice. In 
this system therefore there will be a competition be- 
tween minimization of vortex lattice ordering energy ver- 
sus minimization of the kinetic energy of charged cur- 
rents and magnetic energy. Both the neutronic vortices 
and the rotation-induced vortices in a charged conden- 
sate mixture [15j carry fractional magnetic flux. It is 
known from the other examples 1(1, 2(| that this kind of 
competitions results in formation of composite vortex lat- 
tices if the ordering energy of one of the vortex sublattices 
is low (it is indeed the case in the system in question). 
In particular one neutronic vortex can anchor one or sev- 
eral rotation-induced vortices in the charged condensates 
mixture. There is certainly a range of parameters where 
the resulting composite vortices carry very small mag- 
netic flux. This would result in a different scenario of 
electromagnetic coupling reduction of a rotation-induced 
vortex lattice and magnetic field carrying structures 21 1 . 
Observe also that in a mixture the London law is modi- 
fied dij. 

In conclusion, in the usually assumed picture, a rota- 
tion of a neutronal star induces neutronic vortex lattices, 
and magnetic field produces protonic fluxtubes. It is also 
assumed in the usual picture that no vortices in the pro- 
tonic condensate can be induced by rotation since single- 
component charged condensate reacts to rotation via the 



London law. In this work it was shown that a mixture of 
hadronic superfluids in a neutron star has a qualitatively 
different rotational response if drag strength exceeds a 
certain threshold p pn /Ppp < — 1/2 (microscopic calcula- 
tions suggest that p pn /Ppp can be as large as —0.7) or if 
several charged condensates arc present. In both these 
cases rotation actually does produce phase windings in 
a charged condensate i.e. protonic (or £~-hyperonic) 
vortices. Both these scenarios, may lead to a strong re- 
duction of electromagnetic coupling between a rotation- 
induced vortex lattice and magnetic-flux carrying struc- 
tures. If the resulting interaction strength falls below 
the temperature scale, in some region of a neutron star, 
it can lead to decoupling of superfluid and flux-bearing 
components. 
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